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CHAPTER  ! 


INTRODUCTION 


The  discharge  of  a  compressible  fluid  from  its  pressurized  container 
into  an  environment,  such  as  the  atmosphere,  and  the  inflow  of  the  fluid 
from  the  outside  to  a  vessel,  are  frequently  encountered  in  engineering 
applications.   The  discharge  of  combustion  products  from  their  combustion 
chamber;  the  discharge  of  fuel  from  its  pressurized  tank;  and  flow  processes 
in  refrigerating  and  heating  arc  typ leal  ones . 

The  state  of  the  fluid  inside  a  vessel  may  change  a?  a  result  o.  heat 
flow  and/or  mass  flow  between  the  fluid  and  the  surroundings.   Heat  n.ay  be 
exchanged  between  the  vessel  itself  and  the  fluid  it  contains.   For  mass 
transfer,  the  flow  rate  to  or  from  a  vessel  is  ordinarily  dependent  on  the 
pressure  difference  in  the  duct  and  the  configuration  of  the  duct  that  car- 
ries the  fluid  into  or  out  of  the- vessel ,   A  converging  nozzle  or  a  pipe  is 
usually  used  for  the  discharge  purpose  and  the  kinetic  energy  is  often 
neglected  in  the  analyses  of  inflow  problems.   The  type  of  the  inflow  duet 
is,  therefore,  immaterial. 

The  procedures  as  demonstrated  herein  for  the  discharge  process  provide 
a  way  to  predict  the  conditions  of  a  perfect  gas  inside  a  vessel  from  the 
measurable  elapsed  time.   Also,  the  conditions  can  be  predicted  from  the 
total  amount  of  mass  discharged  if  the  outflow  gas  can  be  weighed  without 
any  difficulties.   As  to  the  inflow  process,  the  way  of  prediction  of  stat 

.','■.'  on   '"      I        or  heat  transferred  and/or  inflow 
gas ,  is  denonstrated - 


CHAPTER    II 

Till:  PROBLEM 

Nomenclature 

English  alphabets: 

A  Cross   sectional    area,    ft 

ft  Integral    constant 

C  Polytropic   constant   defined   by  PV    =   C 

c  Specific   heat  at    constant   pressure,    Btu/lbm   F 

P 

c  Specific   heat  at   constant   volume,    Btu/lbm  F 

v 

F  Frictional    work  per  unit  mass,    Btu/lbm 

g  Gravitational    acceleration,    32.174   ft/sec 

2 

c  Dimensional    constant.    3  2.17^   lbm-f t/    Ibf/sec 

c 

h  Enthalpy  per-unit",   mass,    Btu/lbm 

J  Joule's   unit  converging    factor,    778   ft-ibf/Btu 

i:  c   /c 

P     v 

1  Height  of   a  vessel ,    ft 

K  (n  +   k)/2nk   --  1,    diniensioiiless 

D.fii  Total    mass    inside   a   vessel,    Ibm 

m  mass    flow  rate,    lb;a/sec 

H  Ck   -   1 l/k,    dimensionless 

n  Polytrop Lc    i.     oner 

P  Pressure,    psia. 

Q  Heat,    Btu 


q  Same  as   Q,    except   fcr  a   unit   mass,    Btu/lbm 

o 
R  Gas   constant,    ft-lbm/lbrn   R 

S  Exponent    in    binoir.ial    integral 

2 

S     Allowable  tensile  sUess,  lbf/in 

T     Temperature,   R 

t, ,  t   Thickness  of  a  vessel,  inches 
1    c 

U     Internal  energy  of  a  system,  Btu 

u     Same  as  U,  except  for  a  unit  mass,  Btu/lbm 


w 

s  ,out 
X     (1 


Velocity,  ft/sec        , 

Capacity  of  a  vessel,  ft  ^ 

Volume  of    the  material    of  a   vessel ,    ft 

3 

Specific  volume,    ft    /lbm 


Rotary  shaft  work  out  per  unit   mass,    Btu/lbra 


(P     -  p  ) 

Elevation,    ft 


Greek   letters: 


o< 


2  gc  k 
1 


2      1 

P2 


1) 


k-1 

k 


y 

e 


V/(Ae  ■-<  n  C    ) 


Crai    c^\   Cpt> 


Elapsed    tirre ,    sec 


k+1   1 


2   ,k-l 


y     Kk  <r?i>    j 


^  n 


Subscripts : 

Numerical  subscript:   representing  that  particular  mode 

e   :  Representing  a  gas  condition  at  exits 

b   :  Representing  an  ambient  condition 

*   :  Representing  a  critical  condition 

i   :  Representing  an  initial  condition  inside  a  vessel 

in  :  Representing  an  inflow  gas  condition 

t   :  Representing  a  condition  of  the  material  wall  of  a  vessel 


Statement  of  the  Problem 

The  Discharge  Process.   A  vessel  of  fixed  capacity  with  e   converging 
nozzle  is  initially  filled  with  a  perfect  gas.   Suddenly  the  gas  is  allowed 
to  escape  through  the  nozzle  into  the  atmosphere.   As  a  result,  the  state 
within  the  vessel  changes  continuously  as  long  as  the  inside  gas  is  being 
discharged.   The  effect  of  heat  exchange  between  the  gas  inside  the  vessel 
and  the  surroundings  is  taken  into  account. 

The  relation  between  the  pressure  change  and  the  elapsed  time,  due  to 
mass  flowing  out,  is  first  to  be  determined,  and  a  numerical  example  is  to 
be  shown. 

The  Inflow  Process .   A  theoretical  analysis  on  the  inflow  problem,  in 
which  a  heat  flow-mass  relation  and  a  heat  flow-pressure  relation  are  con- 
sidered, is  demonstrated.   Also  a  numerical  example  for  this  case  is 
presented. 

The  equations  thus  derived  enable  one  to  predict  the  instantaneous 
state  of  a  perfect  gas  inside  a  vessel. 


CHAPTER  111 

ILLUSTRATIVE  PROBLEt-S 

1 .   Fundamental  Governing  Equations  for  Fluid  Flow. 

The  basic  equations  of  compressible  fluid  flow,  known  as  the  energy 
equation  and  the  equation  of  continuity,  can  be  derived  by  considering  two 
arbitrary  coordinates  normal  to  the  flow  stream,  where  the  conservations  of 
both  mass  and  energy  have  to  be  satisfied, 

■  A,  P,  T,  V\    h, 


V  Ai  pi  ji 


ft  \s  _>.  w 

U     /  ^>2  2  2  u  2    '2 


Fixed  Datum 
Fig.  L.   Control  surfaces  for  mass  and  energy  balances. 

Starting  with  the  assumption  of  the  flou  being  steady  and  one  dimen- 
sional, one  can  make  the  mass  balance  between  the  sections  (l)  and  (2)  of 
Fig.  1  from  the  condition  that  the  instantaneous  mass  flow  rate  across  the 
sections  must  be  identical.   Thus 


Vi 


v2 


A  V 


=  ip.  =  Constant    (1) 


or   in   its    differential    form 

**  4   f  -  f     =     0  (la) 

A  V  V 

The   energy    balance   for   fluids    flowing   in   steady  motion,    as   shown  in 
Fig.    1,    is 

q12     *      lU2   "  V    +    (V2   "   Vl)/2acJ    +    (Z2   "  V    g/8cJ    + 

(P     if     -  P     O/j    +   Ws.out/J  (2) 

or   by    the   definition   h=u    +  P   0/J    , 

q12      =      Chj-V    H    (V*   -  V2)/2g.J    +    (Z2-Zl)g/g.J    ♦ 

Ws,out/J  (2a) 

or   in  its   differential    form 

dq     ^      dh   +  VdV/g.J    +  dZ.g/g.J    +  dV.Tssout/J  (2b) 

Equation    (2b)   applies   both   to   frictionai   and   frictionless    flow.      For 
flew  with   friction 

,2 
q12      =       (U2   ~   V    +     J      I>dl'"/J    "  F/J  <3) 

where 

„2 


|     PdD/J      =      (P2V2   -F^l/r   -     | 


irdivj  e.) 


From    (3)    (4),    and    (2),    one  gets 


22  r2 

(V, -Vf)/2g-J +(Z„ -Z,  )g/g-J +Ws,out/J  =  -  VdP/J-F/J  (5) 

2  1  C  2  I  C  vi 


VdV/g-J  *  dZ  g/g-J  +  d  Ws,out/J  +  dr/J  +  IT  dP/J  ■  0         (5a) 

The  equations  thus  derived  are  to  be  applied  to  solve  the  discharge  and 
the  inflow  problems  herein.   Also,  the  equation  of  state,  PV  =  mRT,  is 
introduced  in  solving  the  problems  related  to  a  perfect  gas. 


2 .   Jlnstgaov  Discharge  of  a  Perfect  Gas  from  a  Ves sel  Through  a  Conversing 

Nozzl e. 
' i  >  Reversible  Adiabatic  Flow  of  a  Perfect  Gas  Through  a  Converging  Nozzle. 
Equation  (2b)  can  be  applied  to  this  cast:  and  is  simplified  as 

VdV   /   g  J      =      -  dh  W 

c 

bn 

Substituting  dh  =  cpd'f  =  777TT7  dT  illto   equation    (6)    and   integrating 
betv/een   the  two   sections,    one  obtains 


2g  kRT 


.  4 


(7) 


k-1 
e,    k 


By    the   relations  —  -  =    (— )        and  PIT   =  RT,    and  neglecting  V,    one  gets 

k-1    ,% 


.fa* 


kPtf 


(l-£>"J 


I    k  -   1  l   -   (P    > 


The  mass  flow  rate  m,  from  equation  CI)) 


(7a) 


A  V 
e  e 


(8) 


Exit 

P 
e 

"We 

T 
e 

V 


Discharge  of  fluid  through  a  converging  r.ozxle. 


10 


2  „     k+1 

f  2"    K  „      f     " 

m     =      A 


(ir^i  i!V>k-<f>kJj 


e  1  k 

When  the  flow  at  the  exit  is  in  the  critical  condition  the  pressure 

ratio,  P  /P.  becomes  P  /P.   Here  P  /P  is  called  the  critical  pressure  ratio 
e 

and  is  defined  as  the  ratio  of  the  pressure  at  a  section  of  maximum  flow  rate 

* 

in  the  nozzle  to  the  pressure  at  the  nozzle  inlet.   As  a  result,  P  /P  is 

obtained  by  the  condition  dm/d(P  IV)   -  0,  thus 


P  *  

_e  P_  r_2_-ik-l 

P       ~     P  'k+lJ 


and 


a  v     _     *,,   HE 
IF*     "    A  T  J  v 


T=(S^f/(w'T■k*K, 


CO 


I_    .    (|_,  k      .    _2_  {10) 


Substituting    (9)   and    (11)    into    (7a),    one  obtains 

*  2g    k  }i  2S„k  ^  *   ^ 

V       -     t^-m       -     [^77^       "     tscURT  ] 

The  mass    flow  rate,   m,    is    therefore 


(12) 


(13) 


(14) 


From  ar,  inspection  of  equations  (8a)  and  (13)  one  realizes  that  the 


11 


flow  rate  is  a  function  of  the  inlet  condition  and  the  pressure  ratio  when  k 

is  specified,  and  is  only  a  function  of  the  inlet  condition  if  the  outlet  of 

the  nozzle  is  at  a  critical  condition. 

(i i )   Reversible  Diabatic  Process  Within  a  Vessel  with  Choking  at  the 

Nozzle  Exit  (P.  <  P  =  P*). 
b  —  e 

It   is   permissible    to  assume   that   the   state  at    the   nozzle   inlet  is 
identical   with    that   inside   the  vessel   at  any  point  of   time.      In  addition,    if 
the  process  within   the  vessel    follows  PV     =  Ct    where  n  and  C  are  constants, 
then  one  can  substitute   this    relation   into   equation    (13).      Therefore 


*  n+JL 

1_ 

„2n 


(15) 


dm 


*  "*1 

AjLp2n     d9 


(16) 


But 


and 


V=      (C/P) 


Differentiating  on   both  sides, 


1-n 

1  V  n 

dm       t.      Vd    l-rr)      »     ~  P  dP 

c  VI 

nCn 


(17) 


Since  dm     =     -dm 


12 


therefore 

l-3n 
d0  =  _  __V B   ?r.  dp  (18) 

C*"D**t 


1-n 

6  =  _  ^ p2n  +  t  (18a) 

C2nA*1>U-n) 

The  integration  constant  is  obtained  by  using  the  boundary  condition 
6  =  0  at  P  =  P .  •   Thus 


P.    ** 
[(=*-)  -  1]      for      n  *   1  (19) 


(n-l)A  Tf/Sl! 


and 


e  =       y —  ln  j,  for  n  =  ,  (20) 

The  mass  in  the  vessel  at  any  state  can  be  expressed  in  terms  of  the 
initial  mass,  m  ,  and  the  pressure  ratio  in  the  vessel. 

I 

V        V  ,P  ,n 


V  P.  tf"  i 

(-T^'R 


.  £  )"  (21) 

l   P . 


The  net  mass  discharged  is 


1  -  (|-)n  j  (21a) 


13 


Similarly, 


1/  = 

1 

p.  l- 

(JL,n 

p 

T     = 

T. 

l 

0   " 

(22) 


(23) 


Heat  added  to  the  substance  inside  a  vessel ,  from  the  first  law  of 
thermodynamics ,  is  equal  to  the  sum  of  the  increase  in  internal  energy  of 
the  substance  and  the  energy  flow  out  ace  impanying  the  discharged  mass. 
Hathematically} 

r 

Q.   =   (U  -  U.)  -  |    helm 


I  ji.T.)  -  f   c  Tdm 
v  j  l    J     p   c 


(c  .'fit  -  c  .ti.T.)  -     e  Tdm  (24) 


Substituting  equations    C 17)    (IS)    and    (23)    into    (24),    one  obtains 
c  V       ■  c   V 

Q  =       J-fp.pl.Xfp.p.l 

<in  R      L  iJ        nR      L  iJ 


(P.    -  P)V     c 
i 

R 


L  n  vJ 


(P.    -  P)V 


(P.     -    P)V 
Qin     =     —J for        n  =  1  (25a) 

Q._     =      0  for        n  =  k  (25b) 


14 
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Fig.  3. 


Time  parameter  versus  pressure  ratio  for  sonic  discbarge  of  a 
perfect  gas. 
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3.0- 
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P 


2.0" 
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Fig.  4. 


Specific  volume  ratio  versus  pressure  ratio  for  sonic 
discharge  of  a  perfect  gas. 
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1  . 0 


Fig.  5.   Temperature  ratio  versus  pressure  ratio  for  sonic 
discharge  of  a  perfect  gas. 
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3 . 0 


P. 


2,0 
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6.   Mass  ratio  versus  pressure  ratio  for  sonic  discharge 
of  a  perfect  gas . 
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7-A.      Boat-in   parameter  versus   pressure   ratio    for  sonic 
discharge  of  a  perfect   <^as , 
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Fig.  7-B. 


Heat-out  parameter  versus  pressure  ratio  for 
sonic  discharge  of  a  perfect  gas. 
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Also,    equation    (25)    gives    the   following   information: 
When 


P.  >  p 

1 

heat  is  added 

tf 

n  <  k 

heat  is  out 

if 

n  >  k 

When 


v< 

p 

heat  is  added 

If 

n>  k 

heat  is  out 

if 

n  <  k 

The  equations  thus  derived  are  summarized  in  Fig.  3  to  Fig.  7-B. 

( i i  i  )   Revarsible  Criabatic  Process^, Within  a  Vessel  with  Subsonic  Flow  at 

the  Nozzle  Exit  (P.  =  P  >  P  ). 
b    e 

Assume  the  gas  within  the  vessel  follows  PV   =  C.   Equation  (8a)  there- 
fore becomes 


2g„  k 


2  k^n   k^l.   -  k^l.  ,  h 
k 


^l-f—  *.***[**     -B 


(26) 


C  (k-1) 


m  =   A  o<  P 


2n 


ink  r  p  k  _  ^ 


(27) 


r   2s.  k 


C   (k-1) 


?Jl 


and   (3  =  p 


k-1 
k 
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The  mass   flowing  out,    dm,    during   time  d6   is    then 


k-n        k-1  h 


2»"[pk     _p] 


dm     =     A    c/,  P2n"-  I  P  *     -  0  1       dB  (28) 


Bui 


^  ,        and         lT  =      (-) 


1-n 

dm       =     —r  P   "     dP  (17) 

c  1_ 

n 
n   C 

dm     =      -  dm  (18) 

c 

Equating    (28)    and    (17)    by   the   relation  given  in    (18),   one  gets 

n+k-2nk        k-1  -\ 


60     =     --^P     2nk       lV     -   p] 


dP  (29) 

A  <*  riC" 


n+k-2nk  k-1  -h 

V  J  P     2nk       [  P  k    -  p)       dP  +  c 


(30) 


where 


A  o<  nC  A  o<  nP. 

e  ei 

n   ■''"  k   —   2  nk  k  ~-  1 

Let      M  =  r—         —  ,  N  =  — ; and     S  =  ~k   ,    then  equation    (30)    becomes 

2nk  k 
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e  =  -  r  J  pM  (pN  -  p  )s  dp  +  2  (3D 

The  general  solution  of  the  binomial  integral  in  equation  (31)  is  riot 
available  and  even  a  particular  solution,  corresponding  to  the  given  con- 
stants k  and  n,  is  usually  very  involved.   For  approximation,  the  Simpson 

method  may  be  applied  in  solving  the  integral.   However,  an  exact  solution 

H  -t-  1  M+J 

can  be  obtained  if  and  only  if  the  condition  that  — jr—  +  S  or    -    is 

an  integer  is  satisfied.   For  instance,  if  k  =  1.4,  then  the  condition 

becomes 

J^-1  =  -P-   +  1.25  =   Integer  (32) 

N       4  n 


M  *  l   +   S  =  -r—  *   0.75  =   Integer  (32a) 

N  4  n 

A  number  of  values  of  n  which  will  satisfy  (32)  or  (32a)  are  listed  in 
Tables  1  and  2.   The  Tables  inform  the  possible  region  of  the  exact  solu- 
tions.  Many  solutions  can  be  obtained  in  the  region  I  >  n  ^  -1,  but  no 
exact  solution  can  be  obtained  if  7<  n<  -7.   Some  of  the  exact  solutions 
corresponding  to  the  given  constants  n  and  k  (-1.4)  are  listed  in  Table  3. 
Civ)   Reversible  Adiabatic  Process  Within  a  Vessel  Containing  Air  with 
Subsonic  Flow  at  the  Nozzle  Exit  ( n=k= 1.4). 

A  reversible  diabatic  process,  following  FO"  =  C,  becomes  a  reversible 
adiabatic  process  if  the  exponent  n  is  replaced  by  k.   Therefore,  the 
pressure-time  relation  of  a  reversible  adiabatic  process  inside  a  vessel 
can  be  obtained  from  the  equation  in  Table  3,  where  the  value  of  n  is  1.4. 


TABLE  1 
VALUES    OF  n  FOR   EXACT  SOLUTIONS    OF  EQUATION    (31),    TYPE  1 
k=1.4 
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(H-i-D/N 


1 

2 

3 

4 

5 

6 

7 

8 

9 

HO 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

30 

40 

30 

60 

70 

80 

90 

too 

110 
120 
130 
140 
150 


-1.4000 
-7.0000 
2.3333 
1.0000 
0.6364 
0.4667 
0.3684 
0.3043 
0.2593 
0.2258 
0.2000 
0.1795 
0.1623 
0.1498 
0.1373 
0.1273 
0.1186 
0.1111 
0.1045 
0.0986 
0.0933 
0.0609 
0.0452 
0.0359 
0.0298 
0.0255 
0.0222 
0.0197 
0.0177 
0.0161 
0.0147 
0.0136 
0.0126 
0.0118 


(M+D/N 

0 

-1 

-2 

-3 

-4 

-5 

-6 

-7 

-8 

-9 

-10 

-]  I 

-12 

-13 

-14 

-15 

-16 

-17 

-18 

-19 

-20 

-30 

-40 

-50 

-60 

-70 

-80 

-90 

-100 

-110 

-120 

-130 

-140 

-150 


-1.4000 
-0.7778 
-0.5385 
-0.4118 
-0.3333 
-0.2800 
-0.2414 
-0.2121 
-0.1892 
-0.170? 
-0.1556 
-0.1429 
-0.1321 
-0.1228 
-0.1148 
-0.1077 
-0.1015 
-0.0959 
-0.0909 
-0.0364 
-0.0824 
-0.0424 
-0.0424 
-0.0341 
-0.0286 
-0.0246 
-0.0215 
-0.0192 
-0.0173 
-0.0157 
-0.0144 
-0.0133 
-0.0124 
-0.0116 
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TABLE  2 
VALUES    OF   n   FOR   EXACT  SOLUTIONS    OF    EQUATION    (31),    TYPE   2 


k=1.4 

(M+D/N+S 

r. 

-2.3333 

1 

7.0000 

2 

1.4000 

3 

0.7778 

'• 

0.5385 

5 

0.4118 

6 

0.3333 

7 

0.2800 

8 

0.2414 

9 

0.2121 

10 

0.1892 

11 

0 . ! 707 

12 

0.1556 

13 

0.1429 

14 

0.1321 

15 

0.1228 

16 

0.1148 

17 

0.1077 

18 

0.1015 

19 

.      0.0959 

20 

0.0909 

30 

0.0598 

40 

0.0446 

50 

0.0355 

60 

0.0295 

70 

0.0253 

80 

0.0221 

90 

0.0196 

1  00 

0.0176 

110 

0.0160 

120 

0.0147 

130 

0.0135 

140 

0.0126 

130 

0.0H7 

(M+D/N+S 

n 

0 

-2.3333 

-1 

-1.0000 

-2 

-0.6364 

-3 

-0.4667 

-4 

-0.3684 

-5 

-0.3043 

-6 

-0.2593 

-7 

-0.2258 

-a 

-0.2000 

-9 

-0.1795 

_!0 

-0.1628 

-11 

-0.1489 

-12 

-0.1373 

-13 

-0.1273 

-14 

-0.1186 

■  15 

-0.1111 

-16 

-0.1045 

-17 

-0.0986 

-18 

-0.0933 

-19 

-0.0886 

-20 

-0.0843 

-30 

-0.0569 

-40 

-0.0429 

-50 

-0.0345 

-60 

-0.0288 

-70 

-0.0247 

-80 

-0.0217 

-90 

-0.0193 

-100 

-0.0174 

-110 

-0.0158 

-120 

-0.0145 

-130 

-0.0134 

-140 

-0.0124 

-150 

-0.0116 
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Let 


where 


e    -    X  i3_  _  JL_i  +  2UL  rA  _  JLn 

B  2N   Lp-2N        p-2NJ  4N        Lp-N        p-NJ 

i  i 


3±Il.  r,_lilAl  _  ,„|L^_x 


8N 


Un|r^l-lnIHil] 


-    e.  ♦  &2  ♦  e3 


1  2N   Lp-2N        p-2NJ 

i 


62  4N        Lp-N       p-K-1 


3  p2r    r,    1 1  +  xi  I      ,    I  1  +  XI 

e3  =     sir"  [lrlrnrl  -  Hn^i 


k-l 


Substituting    B    ■  P     k     =  PN     and     X  =   [l   -  f  p"   ]*      into    the  above 
three   equations, .one  obtains 

P      -4N  P     -3N  \  P.    4N  P.    3N  \       2N9 


©<V      =      Up*)  "   <fJ        ]      =   [(E    '        "   (F">      ] 


?2f 


(34) 


P      -2M         P      -N  h  P.    2N  P.    N  %  4N9 

•»a>      ■      [(f>  "(f      ]      -t^J        -03      --^T      (35. 

e  e  J  P     o 
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i+u-ifn       ^[»-'r>V    8N&3 

6(e3)  =  i„ --^  -  m f-^-q  -  —7^  (36) 

i -li -if-)  ]  i-U-cj*-)    ]      3p 


Equations  (34),  (35)  and  (36 )  are  plotted  in  Fig.  8.   Once  the  pressure 
ratio  P  IV   is  specified,  the  corresponding  values  of  the  functions  8^ ,  &2 
and  G„  can  be  found  immediately. 
(v)   Reversible  Isothermal  Process  Within  a  Vessel  Containing  Air  with 

Subsonic  Flow  at  the  Nozzle  Exit  (n=l). 

Following  the  same  procedures  in  the  case  (iv),  one  gets 


P  -N     f     P.N     |     5119 

e(V  =  [(p-e-)   -l]2  =  [(FJ  -^-— -k 


P  -N     |     P.  N     |     3NG 

«v-[(5»)     -i]2  =  [<^>    -i]2-^f7i  (38) 

P   -N     \  P.  N     £     K6 

e(e3)  =  [(f)   -i]2=[(^  -i]2--~V£         (39) 

Equations  (37),  (38)  and  (39)  are  plotted  in  Fig.  9.   The  values  of 

6  ,  ©  and  9,  corresponding  to  a  pressure  ratio  P  /P  can  be  found  providing 

P   P  ,  A  ,  V  and  1/   are  given.   The  total  time  elapsad,  9,  is  simply  the 
i   e   e   i       i 

sum  of  9  ,  6  and  9  . 

Equations  (37)  and  (38)  can  be  neglected  for  an  approximate  solution. 
From  the  relation  9=9+9+9,  one  gets 
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v  R5/2   ,  F.  N 


(39a) 


(vi)   Calculations  on  Heat  Transfer. 

Assume  the  tank  is  perfectly  insulated,  and  the  gas  and  the  walls  of 
the  tank  are  at  the  same  temperature  at  any  instant.   As  the  temperature  of 
the  material  of  the  tank  decreases  from  Tt.  to  T  ,  or  T.  to  T,  heat  is 
transfered  from  the  tank  walls  into  the  gas  within  the  tank  as  a  result  of 
the  decrease  in  the  internal  energy  of  the  walls.   Mathematically, 


d  Q. 


-d  II 


c   m  di- 
pt  t   t 


in   d  T 


(40) 


But 


d  Q-:, 


-h  dm  +  d(mu) 


c  T  dm  +  m  du  +  u  dm 
P 

c  T  dm  +  m  c  dT  +  c  T  dm 
p  v       v 


(41) 


Equating  (40)  and  (41)  yields 


c   m  dT 
pt   t 


(c   -  C  )T  dm  -  m  C'  dT 
p    v  v 


(42) 


By  separating  variables,  equation  (42)  can  be  integrated.  The  integral 
constant  is  determined  subject  to  the  condition  T  =  T.  at  m  =  m. .  The  final 
result  is 


k-1 


1  + 


m,.  c  „ 


1  + 


t   pt  / 


(43) 


Let  £  = 


m   C 
t   pt 


then 


For  a  perfect  gas 


Rearranging  yields 


Combining  (44),  (45)  and  (40),  one  gets 
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m  . c       k-1 


,  i  ♦  (-)£  i 


(43a) 


P_  _i 

P.   m 


(44) 


1  +  C—iz 


1  +  £ 


m  ,  _   k-1 


(45) 


AU   =  C  „  m„   T. 
t      pt   t   i 


' 


1  +  (~)£ 

m. 
l 

i  +  e 


k-1 


(46) 
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Fig.  8.   Functions  of  time  parameter  versus  pressure  ratio 

for  a  subsonic  discharge  of  air,  0.  =0  and  k=n=1.4. 
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1.0 


eq. (39) 


0.5  ' 

0 


I 

0.1 


0.2 


0.3 


0.4 


0.5 


BCSj) 


S(»2) 


6(6.) 


Functions  of  time  parameter  versus  pressure  ratio 
for  a  subsonic  discharge  of  air,  T=constant  and  n=l . 
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3 .   Unsteady  Inflow  of  a  Perfect  Gas  Lo  a  Vessel  with  Constant  Inflow 

Pressure,  Temperature  and  Enthalpy . 

( i )  The  Adiabatic  Process  and  the  Diabatic  Process  within  a  Vessel  at  an 

Initial  Condition  P.  and  T. . 
i jl 

The  processes  are  irreversible  because  of  the  nixing  of  gases  at  dif- 
ferent temperatures.   Considering  the  vessel  as  a  control  volume  and  taking 
an  energy  balance,  one  comes  to  the  result 

Q    =   (U  -  U. )  -  (ra  -  m. )  h. 
xin  1         1   in 


T  [_!5l—  _  i]  _  c  m4    T   [E-  -  1] 
i  Lm.  T.     J    p  i   in  Lm. 


(47) 


Substituting   the  perfect  gas   relation,    PV  =  mRT,    together   with 

c      »  R/J(k  -  1)   and  c      =  kR/J(k  -   1)   into   equation    (47),    one  gets 
v  p 

P.V  „  k  T.         P  T. 

0         =     * [(--  -  1)   -  ~   ( ~  -   1)1  (^8) 

Qin  J(k  -  1)    LtP.  T.         'P.    T        "J 


If  the  process  follows  PV   =  C,  equation  (48)  becomes 

1 


^in  =  J^bt-"--^^'      -1]} 


(49) 


For  an  adiabatic  process  within  a  vessel,  Q,   ■  0,  the  general  equation 
(48)  becomes 


dP 

1 

1 

dT 

P 

T 

T 

k  T. 
in 

(Q    =  0)  (50) 

in 


n                 P  ,T  .n-1 

If  the  process  still  follows  P  V      =  C,  the  relation  —  "  (--)    can  be 

i  i 
expressed  in  its  differential  form  as 
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f   -    ^T   f 


Equating    (50)   and    (51),    one  gets 
n  I 


n  -   1  _       T 

k  T. 
in 


Rearranging  yields 
k  T. 


(52) 


(52a) 


The  equation  (52a)  thus  derived  is  subject  to  two  conditions,  namely, 
0   =0  and  PIT  =  a  constant.   When  the  process  is  isothermal  or  n  *=  1,  the 

\Ltl 

temperature  inside  the  vessel  remains  constant  at  a  value  equal  to  k«T   . 

It  also  ?ives  an  information  that  T  can  not  be  equal  to  T.   because  it  will 
^  n        in 

result  in  n  -  k,  an  irreversible,  diabatic,  isentropic  process,  which  is 

inconsistent  with  the  case  as  discussed  here.   As  shown  in  the  next  section. 

n  =  1  and  T/T.   ■  k  applies  for  the  case  of  a  perfect  gas  flowing  to  an 

initially  evacuated  vessel  in  which  there  is  no  heat  flow. 

( i i )   Adiabatic  Process  Within  a  Vessel  at  an  Initial  Condition  of  Evacuation . 

Since  Q   =0,  U  -  0  (m  =0),  the  energy  equation  (47)  can  be 
in       l       l  oj      -t 

simplified  as 

mc   T.    =  m  c  T  (53) 

pin        v 


-I_  -  Is. 
T.      c 


(53a) 


This  gives  an  interesting  result  that  for  constant  inlet  temperature 
the  states  in  the  vessel  are  isothermal.   Furthermore,  one  can  find  the 
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relation  between  the  mass  flowing  in  and  the  pressure  within  the  vessel  from 
the  perfect  gas  relation.   Thu5 

|j;  =  |  =  Constant  (54) 

( i i i )  Diabatic  Process  Within  a  Vessel  at  an  Initial  Condition  of 
Evacuation. 
Since  U.  =  0  (ra.  =  0),  the  energy  equation  (47)  is  simplified  as 

Q.    =   U  -  mh.  (55) 

in  in 


Qin   "   I? 


r^b-^) 


If  the  process  follows  PI/  =  C,  one  can  substitute  the  relation 

n-1 

P        n 
T  =  T.    (—-)  into   equation    (55a),    thus 


P  V         f ,        ,  Tin 


(56) 


Vi^ 


Equation  (56)  is  applicable  to  the  adiabatic  process  also.   Substitut- 
ing 0   =  0  into  equation  (55a).  one  finds  the  condition  k  =  T/T,   has  to  be 
*  -in  M  in 

satisfied.   Moreover ,  by  comparing  equation  (55a)  with  equation  (56),  one 

finds  T  =--  T.  ,  if  n  =  1.   The  tank  temperature  remains  constant  in  an 
i 

adiabatic  process  if  n  is  unity,  and  if  the  tank  is  initially  evacuated . 


36 


4 .   Numerical  Examples  for  the  Disc har^e  Problem . 

3 
Statement .   A  cylindrical  tank  having  a  volume  of  100  ft  is  initially 

filled  with  air  at  100  psia  and  600  R.   Suddenly  the  air  is  allowed  to 

escape  to  the  atmosphere  (l^t.7  psia)  through  a  frictionluss  converging 

nozzle  of  one  inch  diameter. 

(i )   Tank  Specifications 

The  size  of  the  tank  is  determined  on  the  condition  of  minimum  wall 
material  to  be  used  for  the  fixed  capacity.   Since  the  configuration  of  the 
tank  is  assumed  to  be  cylindrical,  the  size  of  the  tank  can  be  evaluated. 

v  =  it  r2A   »   ioo 

or        I   =      100/Tl  f  2 

A     =      2  7f  r  2   +   2TlVL 
jy.    =     2ft  (2  r  -  100/7T /•   2)     =     0 

Solving   for    )"  ,    one   gets 

I 
P  =      (50/ft  )3      »      2.515    ft. 

and       £=      100/71  T  2     =     5.031    ft. 

If    the  allowable   tensile   stress   of   the  material    being  used   is   18,000 
psi,    the   thickness   of   the  material 

LL.     ..      100  x   2.515   x   12     ~     _1_" 
jt.     '       2S       ~  2  x   18,000  12 

,  pf  100   x   2.515   x    12     w     1" 

and      tc     -     —     -  Ja^oO  =      6 

One  may   select  a   steel    plate   having  a    thickness   of   ^  inch. 
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( i  i  )   Reversible  Adiabatic  Process  Within  the  Tank  ( n  =  k  =  1 . A ) . 

It  is  assumed  that  the  geometry  of  the  converging  nozzle  is  fixed.   The 

characteristic  of  this  kind  of  nozzle  indicates  that  the  state  at  the  nozzle 

exit  is  in  the  critical  condition  (P  =  P  )  whenever  Pfa  <  P  .   On  the  other 

hand,  the  state  at  the  nozzle  exit  becomes  subsonic  with  P  =  Pb  when 

P  ^  P  .   Both  critical  flow  and  subsonic  flow  will  exist  in  this  problem  as 
b  ' 

the  tank  pressure  is  reducing  from  100  psia  to  14.7  psia. 

When  P  =  P  ■  P  =  14.7  psia,  the  corresponding  pressure  in  the  tank 
e    b 


14.7 


[_2_]k-l 
Lk  +  1J 


14.7 
0.52828 


=   27.8261   psia 


Therefore,  the  regions  are  classified  as 

iOO  >  P  .>  27.8261     Critical  condition  at  exit 
27.8261  >  P  >  14.7    Subsonic  condition  at  exit 

For  the  region  100  >  P  >  27.8261,  the  following  equations  apply. 


2V 


(k  -  1)  A  lji/RT. 


P  ** 
[CpMk 


n 


(19) 


Lk  +  1J 


(9) 


A   H>   p2k 

!_ 

-2n 


(15) 


k-1 
T.       P. 


& 


(22) 
(23) 
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For   the  region   27.8261  >  P  ^  14.7,    the   following   equations  apply. 


N  I         „-2N  „-2N  J 


r  r(1  -KN)* 


r"  -Ppi  '      (i  -pp-VI 

I  -N  „-N  J 


4N  p/.  p 


k-1 


■H 


(33) 


~     =      (f->    k  (22) 


P     i 


{.j-r  (23) 


a      ■      A    o<    [P   k      -  p]  (27) 

The  properties  P,  m,  V*  ,  and  T  of  the  air  within  the  tank,  and  the 
elapsed  time  corresponding  to  the  pressure  reduction  are  evaluated  by  the 
above  equations,  and  are  expressed  in  Fig.  10  and  Table  4.   Also,  by  using 
the  charts  ranging  from  Fig.  3  to  Fig.  8,  the  came  answers  can  be  expected. 


TABLE  4 


ANSWERS   TO  THE  SELECTED  DISCHARGE  PROBLEM,    Q. 
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100 

52.8281 

52.8281 

o.oooo 

2.2208 

600.000 

1.7058 

95 

50.1867 

50.1867 

0.9707 

2.3037 

591.271 

1.6324 

90 

47.5453 

47.5453 

2.0016 

2.3944 

582.207 

1.5585 

B5 

44.9039 

44.9039 

3.1002 

2.4942 

572.777 

1.4840 

80 

42.2625 

42.2625 

4.2752 

2.6046 

562.941 

1.4088 

75 

39.6211 

39.6211 

5.5373 

2.7275 

552.655 

1.3330 

70 

36.9797 

36.9797 

6.8995 

2.8652 

541.868 

1.2565 

C3 

34.3383 

34.3383 

8.3777 

3.0210 

530.515 

1.1791 

60 

31.6969 

31.6969 

9.9919 

3.1987 

518.520 

1.1009 

55 

29.0555 

29.0555 

11.7677 

3 . 4039 

505.789 

1.0218 

50 

26.4141  • 

26.4141 

13.7385 

3.6437 

492.201 

0.9412 

45 

23.7727 

23.7727 

15.9484 

3.9285 

477.605 

0.8604 

40 

21.1312 

21.1312 

18.4587 

4.2733 

461.800 

0.7777 

35 

18.4899 

18.4899 

21.3561 

4.7010 

444.514 

0.6939 

30 

15.8484 

15.8484 

24.7704 

5.2481 

425.361 

0.6078 

27.8261 

14.7 

14.7 

26.4641 

5.5378 

416.316 

0.5698 

26 

- 

14.7 

+1.5454 

5.8129 

408.320 

0.5360 

24 

- 

14.7 

+3.4010 

6.1550 

399.088 

0.4941 

22 

- 

14.7 

+5.4849 

6.5496 

389.289 

0.4453 

20 

- 

14.7 

+7.9003 

7.0110 

378.831 

0.3864 

18 

- 

14.7 

+10.8633 

7.5590 

367.597 

0.3111 

16 

- 

14.7 

+15.0377 

8.2225 

355.432 

0.1996 

14.7 

- 

14.7 

+22.1562 

8.7355 

346.930 

0.0000 

40 


4-1  \ 


1                     1 

1                     1 

CO               , 

1 

o 

o 

o 

o 

1 

in 

1 

I          1 

Al 


( i  i  i  )   Reversible  Isothermal  Process  Within  the  Tank  ( n  =  1 ) . 

The  two  pressure  regions  still  exist  in  this  problem  as  the  pressure  is 

reduced  from  P.  to  P  . 
1     e 

For  the  region  100  >  P  >  27.8261,  the  following  equations  apply. 

(P.  -  P)V 

e    =  — — ; (25a) 

m         J 


*  ,  r—  E 

A  -<f/RT. 

•'       l 


P. 
In  -1-  (20) 


P  (15) 


A?, 


Also,  equations  (22),  (23)  and  (9)  apply. 

For  the  region  27.8261  >  P>  14.7,  the  following  equations  apply. 

From  Table  3  (n  =  1 ) 


5 

1  r,PN  -  a 
5N 


Uv!!  -   (J  )z  -  (P«  -  f   )2] 


k-n         I 
&     =  Aec<P2nk  [PN  -  p  ]2  (27) 

Also,  equations  PU  =  C,  and  (25a)  apply. 

The  results  are  summarized  in  Table  5,  and  are  expressed  in  Fig.  11. 
The  charts  ranging  from  Fig.  3  to  Fig.  9  are  available  for  the  calculation 
also. 
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TABLE  5 
ANSWERS  TO  THE  SELECTED  DISCHARGE  PROBLEM 


T=600  R= 

Constant 

P 

* 

1-' 

P 

0 

0 

1/ 

\n 

lii 

100 

52.8281 

52.8281 

0.0000 

2.2208 

0 

1.7058 

9  5 

50.1867 

50.1867 

1.3540 

2.3377 

92.545 

1.6205 

90 

47.5453 

47.5453 

2.7812 

2.4678 

185.090 

1.5352 

85 

44.9039 

44.9039 

4.2901 

2.6127 

277.635 

1.4499 

80 

42.2625 

42.2625 

5.8904 

2.7760 

370.180 

1.3646 

75 

39.6211 

39.6211 

7.5941 

2.9611 

462.725 

1.2793 

7  0 

36.9797 

36.9797 

9.4153 

3.1726 

555.270 

1.1940 

65 

34.3383 

34.3383 

11.3715 

3.4167 

647.815 

1.1088 

60 

31.6969 

31.6969 

13.4845 

3.7014 

740.360 

1.0235 

55 

29.0555 

29.0555 

15.7814 

4.0379 

832.905 

0.9382 

50 

26.4141 

26.4141 

18.2973 

4.4417 

925.450 

0.8529 

45 

23.7727 

23.7727 

2i.0786 

4.9352 

1017.995 

0.7676 

40 

21.1312 

21.1312 

24.1877 

5.5521 

1110.540 

0.6823 

35 

18.4899 

18.4899 

27.7126 

6.3452 

1203.085 

0.5970 

30 

15.8484 

15.8484 

31.7818 

7.4028 

1295.630 

0.5117 

27.8261 

14.7 

14.7 

33.7675 

7.9811 

1335.868 

0.4747 

26 

- 

14.7 

+1.7934 

8.5417 

1368.666 

0.4422 

24 

- 

14.7 

+3.9242 

9.2535 

1406.684 

0 . 4030 

22 

- 

14.7 

+6.2886 

10.0947 

1443.702 

0.3587 

20 

- 

14.7 

+8.9935 

11.1042 

1480.720 

0.3071 

IS 

- 

14.7 

+12.2638 

12.3380 

1517.738 

0.2435 

16 

- 

14.7 

+16.7960 

13.8S02 

1554.756 

0.1536 

14.7 

- 

14.7 

+24.1367 

15.1077 

1578.818 

0.0000 
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( i v )   Reversible  Diabatic,  Non-Isothermal  Process  Within  the  Tank . 

Assume  the  tank  is  perfectly  insulated,  and  the  air  and  tank  tempera- 
tures are  the  same  at  any  instant. 

From  case  (i )  of  this  example,  the  data  of  the  tank  are  5.031  ft.  high, 
2.515  ft.  radius,  and  i,  inch  thickness.   On  these  bases  the  volume  of  the 
material  (steel)  of  the  tank,  V,  is  2.484  cu.  ft.,  and  1213.8  lb.  (the 
density  of  steel  is  assumed  to  be  7.83  x  62.4  =  488.6  lb/ft). 

For  the  steel  walls  of  the  tank,  assume 

c     =   c     =0.12   Btu/lbm°F 
pt      vt 

100  x  144  x  100 

— ^T^ 7^71 x  0.1714 

__53.3  x  600 =    0529S7 

1213.8  x  0.12 

Substituting  £  =  0.053,  and  P/P.  =.  14.7/100  into  equation  (45),  one 
can  solve  for  m/m.  by  trial  and  error.   Thus 

E-  =  0.1496 
m. 
l 

From  equation    (44) 


T  P  mi  14.7  1 


T.  P.    m  100      '    0.1496 

l  i 

Therefore 


0.98262 


j  n 


c    .  rn  [T.    -  T]  =     c    „   m.    T.    [l    -  ~] 

pt      t  l     l             J  pt      L       l    L            1  .  J 

0.12  x  1213.8  x  600  [l  -  0.98262] 

1518.7  Btu. 
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The   answer  is  checked  as   follows. 

U.    =   c   m.T.    =  0.1714  x  '0°',U\'10(1   *  600  =  4630.69      Btu. 
i  v   i    i  53.3   x  600 

U  =  cvmT  =   0.1714  x   [0.1496   x   ^^ 6g010°]   *   1^00  x  0.98262] 

=  680.7      Btu 

h  =  c   Tdra  ■  c   T(m.    -  m) 

out  p  p        1 

„      ,         r60O   +  6O0   x  0. 982,62  n        r100   x   144  x   100       ,,    ,.qfi 
=  0.24   x   [ ■ J ]   *   L      53.3   x  600  "  °A  l96 

100  x   144  x   lOO-i  _  0.24  x   14400  x  0.6504  x  594.78 
x       53.3   x  600        J  319.8 

=   5468.07 

0       »   (U  -  U   )  +  h  =    (680.7    -  4630.69)   +  5468.07   =   1518.1      Btu. 

Un  i  out 
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5 .   Numerical  Examples  for  Inflow  Problems . 

3 
Statement.   A  cylindrical  tank  having  a  volume  of  100  ft  is  initially 

filled  with  air  at  20  psia  and  500  R.   Air  at  100  psia  and  700  R  is  allowed 

to  flow  into  the  tank.   The  inflow  kinetic  energy  is  neglected. 

(i)   Aciiabatic  Process. 

Equation  (48)  applies.   By  the  condition  Q.   =  0,  one  gets 

k  T.    P  T. 

P.  T.    lP.  T 

l  ii 

and  from  the  equation  of  state,  PV  =  mRT,  one  gets 

T     P  m 

—•  = (V  =  Constant) 

I      P  m. 


Combining  the  above  two  equations  yields 


k  T. 
P         in  .m     ,  ,    , 
-  -   ——(--I)  ♦  1 


MJSn7-°°  ("--  1)  »  1   -   1-96  (  — )  -  0.96 
500     m.  ra. 


Rearranging  yields 


»   (f-  +  0.96)     20  x  144  x  100  P..  +  n  o6) 

i 53.3  x  500    20 


and 


1.96  1.96 


PV     P  x  100  x  144     1.96  x  25  P 
mR       ra  x  53.3         P.^^r 

5o~  +  0-96 

The  results  are  shown  in  Fig.  12  and  are  summarized  in  Table  6. 


TABLE  6 
ANSWERS   TO  THE  SELECTED    INFLOW  PROBLEM 
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Q.     »  0, 


10.8068   Ibn. 


Am 


20 

500 

10.8068 

9.2535 

0 

25 

554.298 

12.1852 

8.2057 

1.3784 

30 

597.561 

13.5639 

7.3727 

2.7571 

35 

632.841 

14.9420 

6.6925 

4.1352 

40 

662,162 

16.3204 

6.1273 

5.5136 

45 

686.915 

17.6988 

5.6501 

6.8920 

50 

703.092 

19.0772 

5.2419 

8.2704 

55 

726.415 

20.4556 

4.8886 

10.6488 

60 

742 . 424 

21.8341 

4.5800 

11.0273 

65 

756.532 

23.2125 

4.3080 

12.4057 

70 

769.058 

24.5909 

4.0665 

13.7841 

75 

780.255 

25.9692 

3.8507 

15.1624 

60 

790.322 

27.347  7 

3.6566 

16.5409 

85 

799.424 

28.7261 

3.4812 

17.9193 

90 

807.692 

30.1045 

3.3218 

19.2977 

95 

815.236 

31.4829 

3.1763 

20.6761 

100 

822.148 

32.8614 

3.0431 

22.0546 

48 


100 


600 


700 


800 


I  m  lbm 


-IT  OR 


to 

_1 


ft  /lbm 


Fig.  12.   The  example  of  the  inflow  of  air,  Q.   =  0. 
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(ii)   Pi.  aba  tic  Process  (n  -  constant). 

Assume  n  =  1.25  =  constant.   Therefore  equation  (49)  applies. 

Qm  =  j(k-i)  1(p7  -  n  -  T~"  L(F}    -  1]) 

_i 

20  x  144  x  100  f  .P    .,    1.4  x  700  r.P  .1.25  _  .  -,  1 
778  (1.4-1)  1/20  500    Ll20         JJ 


From  the  general  equation  (47) 


Q       =  in  c      T  -  m.    c     T.    -  rn  c      T,      +  m.    c      T. 
xin  v  ivi  pin  lpin 

=  m    (c      T   •-  c      T.     )    -  m.     (c      T.    -   c      T.     ) 
v  p      in  l        v      i  p      in 


Rearranging  yields 


Q.      +  m.  (c      T.    -  c      T.     ) 
in           l        vi  p      in 

c  T  -  c      T. 
v  pin 


Q.       +  m.     (c      T.    -   c      T.     ) 
ja_n i        vi p      in 

n-1 

c      T.    (E-)   n  -  c      T. 
v      l      P .  pin 


<  144  x  100  |P     ■    1.4  x  700  r  .V    .1.25   . -,  1 
(1.4  -  1 )   I l20  500     L  20      "   J-> 


20 

27  §_ 

0.25 


0.1714  x  500  x  (ffj'1'25  "  °-24  x  700 


^H^-^-  «  (0-1714  *  500  -  0.24  x  700) 
53  .  3  x  500 


The  results  are  shown  in  Fig.  13  and  is  summarized  in  Table  7. 
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n  =   1.2 

5  =  constant 

P 

«in 

20 

0 

.  25 

*-123.144 

30 

-232.284 

3.'. 

-330.207 

40 

-418.821 

45 

-499.509 

50 

-573.325 

53 

-641.087 

60 

-703.452 

65 

-760.957 

70 

-814.046 

75 

-863.097 

80 

-908.430 

t  j 

-950.323 

90 

-989.018 

95 

-1024.727 

Km 

-1057.640 

TABLE  7 
ANSWERS   TO  THE  SELECTED   INFLOW  PROBLEM 


10.8068 
12.9169 
14.9437 
16.9035 
18.8080 
20.6652 
22.4814 
24.2616 
26.0096 
27.7286 
29.4213 
31.0900 
32.7364 
34.3624 
35.9693 
37.5583 
39.1307 


Q.    ' 

in 

'lbm 

0 

-9. 

54 

-15, 

57 

-19. 

53 

-22. 

25 

-24. 

16 

-25 

.50 

-26 

,41 

-27. 

,07 

-27 

,46 

-27 

,62 

-27 

.85 

-27 

,80 

-27 

.75 

-27 

.50 

-27 

.32 

-27 

.00 

Negati ve   sign   indicates    heat-out 
Total   mass    inside   the   tank 
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80 


60 


40 


A 

/ 

' 

§/                         /    I 

/    atf^X            / 

• 

10 

I 


500 


20 


30 
_1_ 


10 


20 


40 


30 

) 


1000 

1 


H  Q  ..  Btu. 


-j  m  lbm 

-t   q   .  Btu/ 
OUC   I  bra 


Fig.  13.   The  example  of  the  inflow  of  air,  n  ■  1.25  =  constant. 
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(iii)      Diabatic   Process    (q.      =   Constant). 

Assume  50   Btu  is  added  per  pound  of  mass    inside   the   tank.      Equation 
(47)   applies. 

Q       =  m   (c      T  -  c      T      )   -  m.    (c      T.    -  C      T.     ) 
vin  V  p      in  l        v      l  p      in 

m . 

q       =    (c     T  -  c      T      )   -  —  (c      T.    -  c      T.     ) 
4in  v  p      in  m  v      l  p     in 

m.RT 

=    (c      T  -  c      T.     )    -  -~-   (c      T.    -  c      T.     ) 
v  p      in  PV  vi  p      in 


Rearranging  yields 


q.      +  c     T. 
in  p     in 


R  m 
cv  ~~vt  (lJv  Ti   "  cp  Tin> 

50  +  0.24  x  700 


0.1714  -  ~l°     _  (0.1714  x  500  -  0.24  x  700) 
500  x  P 


218 


o.nu*5H 


PV     P  x   144  x  100 

RT  "  53  3._ 218 

0.1714  +  ff-f 


A  mean  value  of  n  can  be  evaluated  by  using  the  relation 


n-1 

I-  .   (?-)  " 
T.       F. 


Rearranging  yields 


53 


_! _  l 

,    ,T  .  ,    r 218 T 

In  (— )         In  [ of^       J 

_i_  (0.1714  +  —~)   y.   500 

p       ,  _ ; ii-i 

Fi  ln  L2oJ 


The  results  are  shown  in  Fig.  14  and  are  listed  in  Table  8. 
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TABLE  8 
ANSWERS  TO  THE  SELECTED   INFLOW  PROBLEM 


q.      =   SO   Btu/lbm 


n 


20  648.809  8.3281  0 

25  719.282  9.3902  -1.588 

30  775.433  10.4523  -12.159 

35  821.225  11.5144  8.823 

40  852.283  12.5765  4.570 

45  891.413  13.6386  3.484 

50  918.961  14.7007  2.978 

55  942.584  15.7627  2.180 

60  963.465  16.8243  2.482 

65  981.778  17.8870  2.339 

70  998.038  18.9490  2.231 

75  1018.572  20.0111  2.105 

80  1025.641  21.0732  2.076 

85  1037.456  22.1352  2.018 

90  1048.189  23.1973  i.969 

95  1057.982  24.2594  1.927 

100  1066.954  25.3215  1.890 


55 


500 

750 

1000 

'           1 

0 

10 

20 
1 

0 

1 

? 

10 

1 1 

rii   1  bm 


Fig.    14.      The   example   of    the    inflow  of  air,    q. 
constant. 
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CHAPTER  IV 


CONCLUSION  AND  DISCUSSION 


The  change  of  states  of  a  perfect  gas,  due  to  mass  flow  and/or  heat 
flow,  follows  the  basic  equations  including  PV  =  mRT  and  the  equations 
derived  from  the  first  law  of  thermodynamics . 

In  the  discharge  process  all  the  properties  such  as  mass,  temperature 
and  specific  volume,  and  the  elapsed  time  can  be  expressed  as  a  function 
of  both  the  pressure  ratio,  P./E,  and  the  initial  conditions.   An  ir.stan- 
taneous  state  inside  a  vessel  can,  therefore,  be  predicted  from  the  quantity 
of  the  outgoing  mass  by  utilizing  equation  (21a),  or  can  be  predicted  from 
the  elapsed  time  by  using  equations  (19)  and  (20)  for  a  sonic  discharge 
condition  at  the  nozzle  exit. 

(i)   When  a  perfect  gas  is  discharged  to  an  environment  through  a  con- 
verging nozzle  providing  the  flow  is  in  a  critical  condition,  the  figures 
ranging  from  Fig.  3  to  Fig.  7-B  apply.   When  it  is  discharged  by  means  of  a 
Laval  nozzle,  the  figures  are  available  also,  if  the  flow  at  the  throat  is 
still  in  a  critical  state. 

(ii)   No  general  solution  in  pressure-time  relation  is  available  when 
the  exit  is  in  a  subsonic  discharge  condition.   However,  an  exact  solution 
corresponding  to  a  set  of  n  and  k  is  possible  if  the  condition 
(H  +  1  )/N  +  S  =  an  integer  or  (M  +  1 )/N  =  an  integer  is  satisfied.   Some 
of  the  solutions,  assuming  k  =  1.4,  are  shown  in  Table  3. 

In  tr=  inflow  process,  the  quantity  of  mass  inflow  is  assumed  to  be 
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measurable.   Therefore,  the  prediction  of  the  states  inside  the  tank  can  be 
accomplished  from  the  known  amount  of  mass  flowing  in. 

(iii)  If  a  tank  is  initially  filled  with  gas  having  an  initial  tem- 
perature equal  to  k-T.  ,  the  temperature  inside  the  tank  will  remain  at  a 
constant  value  in  an  adiabatic  process,  when  the  inflow  gas  temperature, 

T   ,  remains  constant, 
in 

(iv)   If  a  tank  is  initially  at  a  condition  of  evacuation,  the  tempera- 
ture inside  the  tank  will  remain  a  constant  value  equal  to  k-T   in  an 
adiabatic  process.   This  is  consistent  with  the  result  described  in  iii. 

The  technique  of  predicting  a  state  within  a  gas  vessel  has  been  shown. 
Charts  relating  the  properties  and  the  elapsed  time  for  various  kinds  of 
gases  at  subsonic  discharge  conditions  can  be  plotted  by  the  equations 
demonstrated  herein,  and  it  is  recommended  that  the  construction  of  the 
graphs  would  be  desirable. 
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A  technique,  based  on  mass  flow  and/or  heat  flow,  for  the  prediction  of 
a  state  of  a  perfect  gas  inside  a  vessel  has  been  formulated  by  utilizing  a 
series  of  fundamental  equations  including  the  equation  of  state,  and  those 
for  mass  and  energy  conservations. 

A  vessel  with  a  converging  nozzle  of  fixed  geometry  has  been  selected 
for  the  study  of  the  discharge  problem.   In  the  analyses  of  inflow  problems, 
the  inlet  stagnation  temperature  was  assumed  to  be  constant. 

In  a  discharge  process  a  general  solution  to  the  equation,  representing 
the  relation  between  the  instantaneous  vessel  pressure  and  the  elapsed  tiiae, 
is  available  for  the  case  of  a  sonic  discharge  at  the  outlet.   Only  a  par- 
ticular solution,  corresponding  to  a  specified  gas  and  heat  flow  process,  is 
available  when  the  discharge  velocity  is  subsonic.   The  analyses,  both  on 
sonic  and  subsonic  discharge,  lead  to  a  result  which  enables  one  to  predict 
the  instantaneous  state  of  the  perfect  gas  inside  a  vessel  from  the  elapsed 
time. 

In  an  inflow  process  the  total  amount  of  mass  flowing  to  a  vessel  is 
assumed  to  be  measurable.   The  instantaneous  properties  are,  therefore, 
expressed  as  a  function  of  mass  inflow.   As  a  result,  a  prediction  of  the 
states  inside  the  gas  vessel  is  possible  from  the  known  quantity  of  mass 
flowing  in. 

Two  numerical  examples  on  the  discharge  and  inflow  of  air  are 
presented. 


